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ABSTRACT 

The paper describes a three-step estimator to calibrate a Three-Axis Magnetometer (TAM) using TAM and 
slit Sun or star sensor measurements. In the first step, the Calibration Utility forms a loss function from the 
residuals of the magnitude of the geomagnetic field. This loss function is minimized with respect to biases, 
scale factors, and nonorthogonality corrections. 

The second step minimizes residuals of the projection of the geomagnetic field onto the spin axis under the 
assumption that spacecraft nutation has been suppressed by a nutation damper. Minimization is done with 
respect to various directions of the body spin axis in the TAM frame. The direction of the spin axis in the 
inertial coordinate system required for the residual computation is assumed to be unchanged with time. It 
is either determined independently using other sensors or included in the estimation parameters. In both 
cases all estimation parameters can be found using simple analytical formulas derived in the paper. 

The last step is to minimize a third loss function formed by residuals of the dot product between the 
geomagnetic field and Sun or star vector with respect to the misalignment angle about the body spin axis. 

The method is illustrated by calibrating TAM for the Fast Auroral Snapshot Explorer (FAST) using in- 
flight TAM and Sun sensor data. The estimated parameters include magnetic biases, scale factors, and 
misalignment angles of the spin axis in the TAM frame. Estimation of the misalignment angle about the 
spin axis was inconclusive since (at least for the selected time interval) the Sun vector was about 15 
degrees from the direction of the spin axis; as a result residuals of the dot product between the geomagnetic 
field and Sun vectors were to a large extent minimized as a by-product of the second step. 

1. INTRODUCTION 

The paper presents a new technique for calibrating a Three- Axis Magnetometer (TAM) on a spin-stabilized 
spacecraft. Contrary to the standard attitude-independent method solving only for magnetic biases, the 
developed 3-Step algorithm also solves for scale factors and misalignments (both orthogonal and non- 
orthogonal). 

Due to the low cost of TAMs there is an increasing interest in using them as attitude sensors for both three- 
axis stabilized and spinning spacecraft. As a rule, ground TAM calibration prior to launch is insufficient, 
and the TAM has to be re-calibrated, using in-flight data. TAM calibration is a well-described technique 1,2 
if three-axis attitude estimates are sufficiently frequent and accurate. However, in the case of spin- 
stabilized spacecrafts these high-frequency attitude estimates are usually not available. As a result, 
attitude-dependent calibration algorithms 1,2 become inapplicable. 

Another TAM calibration procedure conventionally used both for three-axis stabilized spacecraft at an 
early stage of the mission and for spinning spacecraft involves the minimization of magnetic field 
magnitude residuals with respect to magnetometer biases. 3 ' 8 An extension to this approach allowing for 
scale factors and nonothogonality corrections to be solved for has been suggested 9 but orthogonal 
misalignments of the TAM still remained undetermined. 
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Section 2 presents the mathematical details of the three-step calibration algorithm. In Section 3, the 
method is illustrated by using in-flight TAM and Sun sensor data from the Fast Auroral Snapshot Explorer 
(FAST) mission. 


2. THREE-STEP ESTIMATION OF TAM CALIBRATION PARAMETERS 

The three-step TAM calibration algorithm, which (contrary to the standard technique 1,2 ) treats orthogonal 
misalignments separately from scale factors and nonorthogonality corrections, is described in this section. 
The main difference between the two sets of calibration parameters is that errors in the TAM orientation do 
not affect the magnitude of measured magnetic field. For this reason one first has to estimate the deviations 
of scale factors and (if necessary) nonorthogonality corrections from their nominal values and only then 
(after re-adjusting the TAM measurements for the estimated errors) one can calibrate the TAM for the 
orthogonal misalignments. This sequence of calibration steps assures that re-adjustment of TAM 
measurements for orthogonal misalignments does not change magnitude residuals (which have been 
already minimized at the previous step). 

The total (generally nonorthogonal) alignment matrix M is separated into symmetric and orthogonal parts, 
S and O, namely, 


M = OS. (l) 

Diagonal elements of the matrix S account for scale factors whereas off-diagonal elements describe 
nonorthogonality corrections. Note that an arbitrary matrix M can always be decomposed into a product 
O S of symmetric and orthogonal matrices: 


s = 

\Im t M 

(2) 

o = 

MS' 1 . 

(3) 


The three steps of the suggested calibration scheme can be summarized as follows. In Step 1 we estimate 
the symmetric matrix S, together with TAM biases, by minimizing magnitude residuals. Steps 2 and 3 deal 
with estimation of the orthogonal alignment matrix O. Step 2 is used to estimate the direction of the body 
spin axis in the TAM frame by minimizing residuals of the projection of the geomagnetic field onto the 
spin axis, under the assumption that spacecraft nutation has been already suppressed by a nutation damper. 
Finally, Step 3 deals with the misalignment of the TAM frame about the spin axis, which is estimated by 
minimizing residuals of the dot product between the Sun and geomagnetic field vectors. Alternatively one 
could minimize residuals of the angle between the Sun and magnetic field vectors. However, such a choice 
of the loss function has a certain deficiency for orbits with large eccentricity. In fact, to compute the angle, 
one first has to divide the dot product by the magnitude of the geomagnetic field. As a result, errors in the 
computed angle rapidly grow as the spacecraft approaches the orbit apogee. 

The order of Steps 2 and 3 is essential. Since rotations around the spin axis do not affect the projection of 
the geomagnetic field on this axis, re-adjustment of TAM measurements at the last step does not change 
residuals minimized at Step 2. 

The following subsections present the details of each of the three steps. 

Step 1: Minimization of Magnitude Residuals (MMR) 

The raw TAM measurements # raw (t k ) with time tags tk are adjusted via the relation: 
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B TAM (t k ) = s[s raw (t k )-^] , 


(4) 


where b is a vector composed of biases for each of the TAM axes and S is a 3x3 symmetric matrix. If the 
pre-launch calibration provides a nonorthogonal alignment matrix M n0 m> then it can be also factored as 
M n0 m = Onom S n0 m , where O nom is an orthogonal transformation from the nominal TAM frame to the 

body coordinate system (BCS). The symmetric matrix S nom and the pre-launch bias estimate b nom serve 

as a starting point for the matrix S and the vector b . In the first-order approximation, the diagonal 
element Sy of the symmetric matrix S describes a scale factor for the magnetometer axis j, whereas its off- 
diagonal elements Sjp ( j>j'=l,2) are related to angles 4>jj' between the axes j and j' via the relation: 
cos <£jj' - 2 Sjj’ . 

Therefore b and S in Eq. (4) include pre-launch estimates b nom and S n0 m of biases, scale factors, and 
nonorthogonality corrections plus any post-launch adjustments to those nominal values. A general way of 
estimating these adjustments to b nom and S no m is to minimize the loss function 


l=e{ s[s raw (t k )-£] 

k 1 



(5) 


with respect to b and S, using known reference vectors fi 1 ^) in geocentric inertial coordinates (GCI). 


To avoid calibration of magnetic torquer coupling coefficients, only time intervals without noticeable 
spikes in torquer magnetic moment are used. This allows contamination of magnetometer measurements 
due to torquer activity to be treated as noisy corrections to biases. 

After the first step has been accomplished, the nominal (pre-launch) orthogonal alignment matrix O nom is 
used to rotate Eq. (3) into the BCS, 

5 MMR (t k ) = O nom 5 TAM (t k ). (6) 

Note that the vectors #^MR are resolved in the BCS, and superscript ‘MMR’ simply indicates that TAM 

measurements in the TAM frame were obtained via Eq. (4) with the estimated values of b and S 

from the MMR. (Since calibration for misalignments of the TAM relative to the body axes has not been 
done yet, it is assumed that the TAM frame in question coincides with the nominal TAM frame defined by 
the alignment matrix O n0 m •) TAM measurements adjusted in this way serve as the starting point for Steps 
2 and 3 to determine the TAM orientation O m i s relative to the nominal TAM frame. 


Step 2: Estimation of Spin-Axis Misalignment in TAM Frame 

The final adjustment of TAM observations relative to the BCS is done via the relation: 

fl B (t k ) = Ofl TAM (t k ), (7) 
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where O — O n om Omis • 

Let us represent the orthogonal matrix O as 

( 8 ) 

4 TAM in the TAM frame, 

The second loss function is introduced as 

L'(Ab,s 3 ,4 TAM ;n s I pin ) = Z [s 3 5 TAM (t k ).P z TAM -Ab-fi 1 ^).^] 2 (9) 

k 

which is minimized with respect to the scalar bias correction Ab, the scale factor S 3 , and the direction 
P™ of the spin axis in the TAM frame, under the assumption that the direction of the spin axis in GCI , 

^spin ’ remains unchanged with time and has been independently estimated either using another pair of 

sensors or by combining partially calibrated TAM measurement vectors in Eq. ( 6 ) with Sun (or star) sensor 

measurements. Alternatively, one can treat the projections of the vector ^spin as additional estimation 

parameters and, in this case, Step 2 combines TAM-only single-axis attitude estimation with a partial TAM 
calibration. 

Including the scale factor S 3 into the set of calibration parameters estimated in Step 2 makes it possible to 

solve the problem in an analytical form, whether «spin * s a ^ nown or solved-for vector . The appropriate 

analytical formulas are discussed in the Appendix. Assuming that Step 1 is sufficiently accurate, then both 
the scale factor correction AS 3 = S 3- 1 and the bias correction Ab are expected to be negligibly small. If 
this is not the case, then one has to re-investigate the output from Step 1 and adjust the third diagonal 

element of the matrix S and the z-component of the vector S b for the corrections AS 3 and Ab from Step 2 . 
One can then minimize the loss function in Eq. (5) only with respect to the two other scale factors and the 

x- and y-projections of the vector S b , and repeat Step 2 . If necessary, additional iterations could be 
performed. 

Note that the loss function in Eq. (9) explicitly assumes that the GCI direction of the body Z axis remains 
unchanged with time. This implies that the spacecraft is in the nominal mode and thus undergoes pure 
rotation. An extension of this algorithm to nonzero nutation angles is discussed in Ref. 10. 

An intermediate set of partially adjusted vector measurements in the Spin Axis Frame (SAF) 

B SAF (t k ) = A 123 ((M,0)B TAM (t k ) 00) 

is then introduced by defining the misalignment angles (j) and 9 via the relation: 

P™ = [sin0, -cos0sin<{>, cos0cos<|)] r . (11) 

The matrix Ai23(^>9>0) in Eq. (10) is the 1-2-3 sequence of Euler rotations with angles (J), 0, and 0. The 
full transformation matrix o from the TAM frame to BCS can be written: 


0 = [P? AM ,Pj AM ,P? AM ]. 

The purpose of Step 2 is to locate a misaligned direction of the body spin axis 
while Step 3 estimates the misalignment component about this direction. 
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O = Ai23(0,0,V) =A3(y) Ai23(^0,O) . 


( 12 ) 


As the zeroth-order approximation to the third Euler angle \J /, one can use its nominal value \|/ n0 m , found 
from the 1-2-3 decomposition Ai 23 ((|>nom >0 nom ^ nom) of the matrix O n om- The post-launch calibration 
of this misalignment component is discussed in Step 3. 


Step 3: Estimation of Misalignment of TAM Frame About Spin Axis 

The last step, which deals with residuals of the dot product between the Sun (or star) and magnetic field 
vectors, requires knowledge of accurate time-tags of the Sun (or star) slit sensor measurements. The 
pertinent loss function has the form: 

L' = E [S B (t>A , (13) 

£ 


where S B (t^) is a Sun measurement vector in BCS, with time tags , 5 T (t^) is the reference Sun 

vector in GCI, B (tg) and B (t^) are magnetic field reference and partially adjusted (via Eqs. (4) 
and (10) above) TAM measurement vectors interpolated to the times of the Sun sensor measurements. 
Minimization of the loss function in Eq. (13) with respect to the Euler angle V|/ is performed numerically. 

Note that Step 3 fails when the Sun vector is either parallel or anti-parallel to the spin axis. One should 
thus expect large errors in estimated values of the Euler angle \j / when the Sun vector becomes 
perpendicular to the body XY -plane. 


3. TAM CALIBRATION USING FAST IN-FLIGHT DATA 

In this section we describe the results of applying the three-step calibration algorithm to actual in-flight 
data from the FAST spacecraft. The FAST is a spin-stabilized spacecraft with a spin rate of 12 revolutions 
per minute (RPM) about the Z-axis. It has an eccentric 350x4200 kilometers orbit around the Earth, with an 
83-degrees (deg) inclination. The nominal direction of the spin axis is 3.5 deg south of the negative orbit 
normal. The onboard sensors include one Horizon Crossing Indicator (HCI), one Spinning Sun Sensor 
(SSS) of 0.5 deg resolution, one TAM of 1 milligauss (mG) resolution, and two electromagnetic torquers. 

We used the 4-orbit time span of FAST data, described in more detail in a previous paper 11 on TAM-only 
attitude/rate determination for spin-stabilized spacecrafts. Since an analysis of torquer magnetic moments 
revealed no pronounced spikes, we disregarded any contamination of TAM measurements due to torquer 
activity. 

The TAM calibration was done in the order specified in the previous section. First, the TAM was 
calibrated for biases and scale factors by minimizing the loss function in Eq. (5). (Pre-launch estimates are 

[0 0 0] for b nom and the identity matrix for M nom .) As a necessary condition for this step to be 

successful, we required that the resultant magnitude residuals form a band with a nearly zero mean and a 
width of about 5 mG. The first step was also supposed to suppress pronounced orbit-dependent variations 
of magnitude residuals, which were observed for magnitude residuals of uncalibrated TAM measurements 
depicted in Fig. la. In addition, since the body Z-bias and the scale factor for the Z-component of TAM 
measurements were also estimated in Step 2, the calibration can be considered sufficiently accurate only if 
the corrections to the bias and the scale factor from Step 2 are small. 
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Figure 1. Magnitude Residuals of Measured Geomagnetic Field 

One should also check the sensitivity of the Step 2 calibration parameters to the choice of the single-axis 

attitude solution for the direction of the spin axis in GCI , ftspin * We use< ^ single-axis attitude solutions 

from the two conventional single-axis attitude estimation schemes: the SSS/HCI Earth-In Method 12 and the 
SSS/TAM Cone Method, 13 as well as instantaneous directions of the spin axis associated with three-axis 
attitude solutions from the SSS/TAM TRIAD. 14 Since both the SSS/TAM Cone Method and the SSS/TAM 
TRIAD initially use adjusted TAM measurements from Step 1, one has to assess how the Step 2 
calibration affects the accuracy of these solutions. 

Finally, one has to analyze residuals of the dot product between the Sun and geomagnetic field vectors. For 
this step it is essential that the SSS provides both the Sun angle (i.e., the angle between the Sun vector and 
the spacecraft Z-axis) and the time when the Sun entered the field of view (FOV) of the SSS. Using these 
observations, one can construct the necessary Sun vector measurements. 

Evaluation of the three-step calibration algorithm was done for the TAM adjustment options listed in 
Table 1. First, we minimized the loss function in Eq. (5) with respect to only biases, since magnitude 
residuals are usually not sensitive to errors in scale factors and nonorthogonality corrections. However, as 
seen from Figures la and lb, minimization of the Step 1 loss function only with respect to biases made 
errors in magnitude residuals near the orbit perigee even more pronounced, while giving an anomalously 
large value of -29.5 mG for b z (Case b in Table 1). Magnitude residuals were uniformly reduced (see 
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Figure lc) only after the loss function was minimized with respect to both biases and scale factors (Case c 
in Table 1) . An additional minimization of the Step 1 loss function with respect to off-diagonal elements of 
the matrix S did not produce reliable results and will be omitted from this discussion. 


Table 1. TAM Adjustment Options 


Description 

Case 


[Sn, S22, S33] 

■a 


no adjustment 

a 

[0, 0 , 0] 

[i, 1, i] 

0, 

0 

only for biases 
from Step 1 

m 

[-0.47, 1.62,-29.52] 

[l.i.i] 

0 , 

0 

for biases and scale 
factors from Step 1 

H 

[-0.74,-2.01, 5.77] 

[1.041, 1.022, 1.024] 

0, 

0 

corrections to Case c 
from Step 2 (TRIAD) 

m 

[-0.74, -2.01, 4.27] 

[1.041, 1.022, 1.039] 

0, 

0 


e 

[-0.74,-2.01, 4.34] 

[1.041, 1.022, 1.032] 

0.70, 

-0.46 


Cases d and e reflect re-adjustments based on calibration results from Step 2, which are summarized in 
Table 2. Since the bias and scale factor corrections Ab and AS3 turned out to be relatively large, we 
adjusted both the bias b z and the diagonal element S33 for these corrections, re-adjusted TAM 
measurements via Eqs. (4) and (6), repeated Step 2, and confirmed that this re-adjustment did not affect 
previously estimated values of misalignment angles (j) and 0. 

The algorithm was implemented for the following two approaches described in the Appendix: 

i) the GCI direction of the spin axis, , is independently estimated using other sensors; 

ii) the GCI direction of the spin axis, nlp[ n , is estimated together with (j) and 0, by minimizing the loss 
function in Eq. (10). 


Table 2. Step 2 Calibration Results Using Different Attitude Estimation Methods 


TAM Adjustment 
Option 

Source of Spin Axis 

Ab (mG) 

As 3 

<J>, 0 (deg) 

Case c 

SSS/TAM (Cone Method) 



0.69, -0.46 

Case c 

SSS/TAM (TRIAD) 

-1.56 

■KKTQHi 

0.69, -0.46 

Case c 

TAM-only 

-1.48 

0.013 

0.67, -0.28 

Case d 

SSS/TAM (TRIAD) 



0.70, -0.46 

Case d 

TAM-only 



0.70, -0.46 

Case e 

TAM-only 

0.00 


0.70, -0.46 


Implementation of the first option was done using the three attitude estimation methods mentioned above: 
the SSS/HCI Earth-In Method, 12 the SSS/TAM Cone Method, 13 and the SSS/TAM TRIAD. 14 Case d in 

Table 1 describes the bias and scale factor re-adjustments based on the bias and scale factor corrections Ab 

and AS3 from Step 2, which was implemented using the TRIAD solution for ^spin » w Rh Case c 

option for adjusted TAM measurements. Case e corresponds to the final output from the TAM-only 
version of the algorithm, which practically converged after the first iteration. 

Note that both the TRIAD and the Step 3 calibration require interpolation of TAM measurement vectors to 
times of Sun sensor measurements. Since FAST is doing 12 revolutions per minute, X and Y components 
of the measured geomagnetic field rapidly oscillate, and the standard interpolation schemes do not work. 
For FAST we took advantage of the fact that the spacecraft rotates with a nearly constant rate, by 
interpolating projections of measured vectors onto the axes of the so-called ‘despun’ frame (which rotates 
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around the Z axis with the frequency of 12 Hertz) . The interpolated vectors were then rotated back to the 
BCS. 

Table 3 compares mean GCI directions of the spin axis estimated using different attitude estimation 
methods. While predictions of the two single-axis attitude estimation methods for the GCI direction of the 
spin axis are relatively close, the TRIAD solution points about half a degree away. One might assume that 
discrepancies come from a TAM misalignment around the spin axis, which would result in large attitude 
errors re-distributed among all three axes. However, as discussed below, this is not the case, and TAM 
calibration has a very little effect on the estimated direction of the spin axis. 


Table 3. Mean GCI Directions of Spin Axis From Different Attitude Estimation Methods 


mam 

Attitude Estimation Method 

Right ascension (deg) 

Declination (deg) 

N/A 

SSS/HCI 

12.79 

-11.34 

Case c 

SSS/TAM (Cone Method) 

12.66 

-11.24 

Cases c and d 

SSS/TAM (TRIAD) 

13.03 

-10.86 

Case c 

TAM-only 

13.16 

-10.71 

Case d 

TAM-only 

13.12 

-10.75 

Case e 

SSS/TAM (TRIAD) 

13.05 

-10.86 

Case e 

TAM-only 

13.16 

-10.75 


Figure 2 depicts instantaneous GCI directions of the spin axis for each of the three methods, using adjusted 
TAM measurements from Case c. It is interesting that the line formed by solutions from the SS S/TAM 
Cone Model (despite a large spread of points along the line) goes directly through the spot formed by 
solutions from the SSS/HCI Earth-In Cone Method, whereas the spot formed by TRIAD solutions lies 
noticeably apart. As discussed below, calibration of TAM measurements for orthogonal misalignments 
reduced deviations of instantaneous directions of the spin axis from its mean position but had no significant 
effect on the estimated mean position. 

The two attitude estimation schemes, the SSS/HCI Earth-In Method and the SSS/TAM TRIAD, also 
predict completely different patterns for deviations of the instantaneous direction of the spin axis from its 
mean position. As seen from Figure 3, the SSS/HCI Earth-In Method predicts a regular orbit-correlated 
motion of the spin axis around its mean position with relatively small nutation, whereas points associated 
with TRIAD solutions form a band. Step 2 calibration (using the Case e adjustment option) slightly 
reduces the band width, but does not change the pattern. 

The simplest explanation for the discrepancy between the mean GCI directions of the spin axis predicted by 
the Single- Axis Attitude Estimation Methods and the TRIAD method is that a measured Sun angle has an 
undetected bias. This bias would have a much smaller effect on TRIAD solutions, because they are based 
on all three components of the measured Sun vector. Such an explanation is supported by the fact (see 
Table 3) that the T AM-only solution for the mean GCI direction of the spin axis is very close to that 
predicted by TRIAD. Bearing in mind that the SSS/HCI Earth-In Method 12 is an integral part of the 
application 15 used in Flight Dynamics operations for single-axis attitude validation, the discrepancy 
requires a more thorough investigation. 

Since the calibration results given by Case e are based only on TAM measurements, this adjustment option 
was selected as a starting point for Step 3. Figure 4 illustrates the effect of calibration on residuals of the 
angle between geomagnetic field and the Sun vector. It is remarkable that these residuals became 
drastically smaller after Step 2, despite the fact that the loss function in Eq. (13) had not yet been 
minimized. Since Sun measurements have not been used for Case e, the small angle residuals serve as an 
additional argument in support of the calibration results. 
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Figure 2. Estimated GCI directions of the spin axis using three different models: 

a) SSS/HCI Earth-In Method, b) SSS/TAM Cone Method, c) SSS/TAM TRIAD 
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Figure 3. Angle between ‘instantaneous’ direction of the spin axis and its mean position 
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Figure 4. Residuals of the angle between the Sun vector and geomagnetic field 

An attempt to further minimize these residuals using the Step 3 loss function did not produce any 
significant improvements in the resultant plot for the residuals. One should be careful, however, with final 
conclusions concerning the magnitude of the misalignment angle, taking into account that, during the 
selected time period, the Sun vector forms a relatively small angle (about 15 degrees) with the spin axis. 

As a result, the main contributions to the dot product between the geomagnetic field and the Sun vector 
come from the projection of the geomagnetic field on the spin axis, whose residuals were indeed minimized 
in Step 2. In other words, the relatively small angle between the Sun vector and the spin axis for this 
FAST data set negatively affects sensitivity of the Step 3 calibration. To clarify this point, Figure 5 depicts 
the RMS of the dot product between Sun and geomagnetic field vectors. 


RMS = 7 L '(V)/n (14) 

as a function of all possible values of the misalignment angle y, where n is the number of good Sun 
measurements. As seen from the figure, there is a shallow minimum near \\f ~ 0, which is associated with 
uncertainties of a few tenths of degree in estimated angles of \|/. This implies that, for this FAST data set. 
Step 3 can estimate the TAM misalignment component about the spin axis only with the accuracy of about 
±0.5 degree. Case e can be thus considered as the final set of calibration parameters, which satisfies all of 
the specified validation criteria (except for the fact that TRIAD and TAM-only predictions for the mean 
GCI direction of the spin axis do not match those from Single-Axis Attitude Estimation Methods.) 


4. CONCLUSIONS 

The main purpose of the paper is to present a new attitude-independent TAM calibration algorithm for 
spin-stabilized spacecraft in the nominal mode and to evaluate its effectiveness, using spacecraft in-flight 
data. The most important result of the paper is a new tool for the attitude estimation of TAM orthogonal 
misalignments, which is achieved in Steps 2 and 3 of the developed 3-Step algorithm. Another important 
development is a Step 1 extension of the conventional scheme for attitude-independent estimation of 
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Figure 5. RMS of residuals of the dot-product between the Sun vector and geomagnetic field. 

magnetic biases to include scale factors and nonorthogonality corrections into a set of estimated calibration 
parameters. 

The algorithm was evaluated, using FAST in-flight data by comparing magnitudes of pre- and post- 
calibration TAM-only measurement residuals. Additional evaluation of the algorithm consisted of 
analyzing mutual consistency of the results from different calibration scenarios, as well as examining an 
effect of the TAM-only calibration on residuals of the angle between the measured Sun and geomagnetic 
field vectors . 

Using these evaluation schemes, we concluded that the FAST TAM was successfully calibrated for three 
biases, three scale factor components, and two components of the orthogonal misalignment. It was 
demonstrated that the calibration only for biases exaggerated the relatively large errors in magnitude 
residuals near the orbit perigee and gave an anomalously large value for the Z-bias. This problem was 
eliminated after minimizing magnitude residuals with respect to both biases and scale factors. (Usually 
Step 1 allows one to estimate only biases. We found that the method was still not sensitive enough for 
reliable estimation of nonorthogonality corrections and the latter were simply neglected.) Based on results 
from Step 2, we also had to improve the body Z-bias and the scale factor for the third TAM axis. It was 
verified that these corrections did not have any noticeable effect on plotted magnitude residuals. 

We found that the output from Step 2 (after the cited re-adjustment of the Z-component of TAM 
measurements) automatically took care of residuals of the angle between the Sun and magnetic field 
vectors. Solving for the misalignment about the spin axis in Step 3 thus turned out to be inconclusive 
(based on our evaluation data set). As explained in the paper, one of the reasons for such a relatively low 
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sensitivity of Step 3 is that the Sun vector is just 15 degrees from the spin axis during the selected time 
period. 

We conclude that the new algorithm provides an effective tool for the in-flight calibration of TAMs on 
spin-stabilized spacecrafts in the nominal mode. It is recommended that the suggested three-step 
algorithm be used when TAM calibration is desired. 

One remaining concern is the observed discrepancy between the mean GCI direction of the spin axis 
estimated via the Single Axis Attitude Estimation Methods and that estimated via TRIAD. This unresolved 
issue deserves more study. 

Other suggestions for further study include: 

(1) the extension 10 of the Step 2 algorithm making it possible to include nutation of an axially-symmetric 
spacecraft into the Step 2 loss function; (2) study of the sensitivity of Step 3 results at larger values of the 
angle between the Sun and magnetic field vectors. 
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APPENDIX 

This appendix provides analytical formulas to minimize the loss function in Eq. (9) with respect to a bias 
correction Ab, the scale factor S 3 , and the directions /^AM ^ body z axis m the TAM frame for the 

known direction /^ pin of the spin axis in the inertial coordinate system. If necessary, one can then 

substitute the derived explicit expressions of Ab, S 3 , and P Z TAM in terms of n* pin into the Step 2 loss 

function and find (also analytically) the absolute minimum of the latter function with respect to n* in . 


Let us first exclude the bias correction Ab from Step 2 loss function, by nullifying the first derivative of V 
with respect to Ab and then using the resultant relation to explicitly express Ab in terms of other 
estimation parameters: 


Ab = mean{ B TAM (t k ) } • P Z TAM - mean{ B 1 (t k ) } • nj pin . 


(Al) 


The new loss function thus takes the form: 


r / / c AT AM , -jl \ y r r a zjTAM \ _ AT AM * nl \ m I "| 2 

L ( S 3>P Z >^spin) — 2 L s 3^^ (^k ) # (^k) # ^spinJ 9 


where we put 

AB tam (t k ) s B TAM (t k ) - mean{ fi TAM (t k ) } 


and 


AB 1 (t k ) s B 1 (t k ) - mean t B l (t k ) } . 


(A2) 


(A3) 

(A4) 
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Let us also introduce three auxiliary 3x3 submatrices: 

{f-TT }jj' = I AbJ am (t k ) AB} am (t k ) , 
k 

{fnlff = LABT AM (t k )AB}(t k ), 

k 

and 

{fn )jj' = £ ABj(t k )ABj'(t k ), 

k 

and represent the loss function given in Eq. (A4) in a simplified form: 

t '( pTAM -I \ _ pTAM f pTAM , -I # f 9 pTAM n -I 

L( T Z ’^spin^-^z #1 TT / z +n spin #1 n n spin ' 2/ z M Tl' z spin' 


(A5) 

(A6) 

(A7) 

(A8) 


For each unit vector n* pin , loss function (A8) has only one extremum with respect to given by the 

linear equation: 


ST AM 


. z - f-^x f TI«spii 


(A9) 


If the direction ^gpj n is known, Eq. (A9) gives the estimation parameters in question: = | />TAM | an( j 

3’ 


pTAM = pTAM 


To minimize the loss function also with respect to n] pin , one simply needs to substitute Eq. (A9) into 
loss function (A8), which gives 

l'(4o) = 4o*f4V (A10) 

where we put 


F=f n -fxifTT f TI 


(All) 


The direction fl^pin sought for must be either parallel or anti-parallel to the eigenvector n s * Din of the matrix 


'spin 


F associated with the lowest eigenvalue; that is, 

I . 


^spin “ - ^spin ^ * ^spin 1 


(A 12) 


The sign in the latter expression should be chosen in such a way that the vector ,P Z TAM found via 
Eq. (A12) forms an acute angle with the nominal direction of the body axis Z in the TAM frame (given by 
the third row of the nominal alignment matrix O n om)- Finally the bias correction Ab is given by Eq. 
(Al). 
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